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Abstract 

Background: Illegal hunting is one of the major threats to vertebrate populations in tropical regions. This 
unsustainable practice has serious consequences not only for the target populations, but also for the dynamics and 
structure of tropical ecosystems. Generally, in cases of suspected illegal hunting, the only evidence available is 
pieces of meat, skin or bone. In these cases, species identification can only be reliably determined using molecular 
technologies. Here, we reported an investigative study of three cases of suspected wildlife poaching in which 
molecular biology techniques were employed to identify the hunted species from remains of meat. 

Findings: By applying cytochrome b (cyt-b) and cytochrome oxidase subunit I (COI) molecular markers, the 
suspected illegal poaching was confirmed by the identification of three wild species, capybara {Hydrochoerus 
hydrochoeris), Chaco Chachalaca {Ortalis conicollis) and Pampas deer {Ozotoceros bezoorticus). In Brazil, hunting is a 
criminal offense, and based on this evidence, the defendants were found guilty and punished with fines; they may 
still be sentenced to prison for a period of 6 to 12 months. 

Conclusions: The genetic analysis used in this investigative study was suitable to diagnose the species killed and 
solve these criminal investigations. Molecular forensic techniques can therefore provide an important tool that 
enables local law enforcement agencies to apprehend illegal poachers. 
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Findings 

It is estimated that illegal hunting kills millions of verte- 
brates per year in tropical rainforests [1]. A number of 
reports have shown that the volume of wild game har- 
vested is unsustainable and has led to the local extinc- 
tion of several populations [1,2]. The use of wild foods 
or the bush meat crisis is one of the major challenges 
for the conservation of large-bodied vertebrates. Poor 
local law enforcement and corruption allow the hunting 
of large vertebrates to continue to be widespread even in 
protected areas [2]. 

Even when hunters are captured, the precise identifica- 
tion of bush meat is often questionable. Moreover, the 
remains of meat, fur, skin and bone are the only evi- 
dence recovered from the crime scene. In such cases, 
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species identification can only be resolved with molecu- 
lar tools [3-5]. 

Here, we reported an investigative study of three sus- 
pected offenses of wildlife poaching. In July 2010, based 
on suspected illegal hunting, a wildlife inspector of the 
Brazilian Environmental Agency (IBAMA) seized and 
sent us biltong samples of a mammal species (MAM1) 
(Figure 1). According to the inspector, the suspect 
claimed that the meat was pork, but there was no evi- 
dence to confirm this assertion (case 1). Along with 
MAM1, we received a second meat sample, which was 
removed from the wings of an unidentified bird species 
(BIRD; case 2) (Figure 1). In February 2011, we received 
a mammal meat sample (MAM2) taken by another wild- 
life inspector of IBAMA. The seized meat was confis- 
cated from the suspects freezer during routine 
surveillance activity (case 3). All three seizures were per- 
formed in the central-western region of Brazil, and spe- 
cies identifications were not possible from morphological 
data. Therefore, the cases could not be characterized as 
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environmental crimes. To enable the appropriate law en- 
forcement actions to be taken, we used the most widely 
used DNA markers for species identification, cyt-Z? [6] 
and COI [7]. DNA was obtained using a standard 
method of phenol-chloroform extraction and precipita- 
tion with ethanol [8]. A negative control was included 
in each DNA extraction process. A portion of about 
1000 bp of cyt-b was amplified using the universal 
primers L14841 (5AAAAAGCTTCCATCCAACATCT- 
CAGCATGAAA3') [9] and H15915 (5AACTGCAGT- 
CATCTCCGGTTTACAAGAC3') [10]. The COI fragment 
(-700 bp) was amplified using the universal primers 
LCO1490 (5' GGTCAACAAATCATAAAGATATTGG3') 
and HC02198 (5' TAAACTTCAGGGTGACCAAA- 
AAATCA3') [11]. The amplifications of both genes were 
carried out in a final volume of 30 ul containing 100 ng of 
DNA, 0.3 mM dNTPs, lx PCR buffer (20 mM Tris-HCl 
pH 8.4, 50 mM KC1), 2.5 mM MgCl 2 , 2.5 U Taq DNA 
polymerase (Invitrogen) and 8 pmol of each primer. The 
PCR conditions were initiated by a denaturation at 94°C 
for 5 min, followed by 40 cycles at 94°C/30 s, 50°C/45 s, 
72°C/45 s and a final extension at 72°C for 10 min. PCR 
contamination controls were performed. Each sample was 
sequenced twice on both strands using the Big Dye v3.1 
Terminator Cycle Sequencing Kit (Applied Biosystems) in 
an ABI 3730 DNA Analyzer (Applied Biosystems). 

The sequences were manually edited using the soft- 
ware BioEdit Sequence Alignment Editor 7.0.9.0 [12]. 
This software was also used to align both strands and 
obtain consensus sequences for each sample. The 
sequences obtained were submitted as independent en- 
tries in a BLAST search for the most similar sequences 
using the default Megablast algorithm parameters [13]. 
For the species diagnosis, we considered the percentage 
homology between query and reference sequence pairs. 
In order to minimize the chance of incorrect species as- 
signment, we adopted a threshold identity value of 
>98 % between the sequences [4]. Since this study aimed 
at the taxonomic identification of the seized meat 



samples, without any prior evidence or indication of the 
species, we regarded uploading the sequences to 
GenBank inadequate. The query sequences can be 
obtained directly from the authors. 

The sequences of the seized samples were compared 
against those of species that were likely to be hunted or 
consumed at the seizure region (Table 1), which were 
downloaded from GenBank. When available, we used 
reference sequences of at least five different specimens 
of each species. The reference and suspected sequences 
were aligned using Clustal W implemented in BioEdit. 
Pairwise genetic distances were obtained according to 
the Kimura two-parameter model [14], and neighbor- 
joining analysis [15] was performed using Mega 3.1 [16]. 
Bootstrap values of the branch configuration in the trees 
were estimated using 1,000 replicates. 

The cyt-b sequencing of BIRD and mammal species 
(MAM1 and MAM2) produced readable sequences of 
approximately 550 and 1100 bp, respectively, which did 
not present insertions, deletions or stop codons. BLAST 
analysis of cyt-b indicated that the MAM1 meat sample 
was highly similar (99%) to the capybara (Hydrochoerus 
hydrochaeris) reference sequence, the MAM2 was simi- 
lar (98%) to Pampas deer (Ozotoceros bezoarticus), and 
the BIRD sample was similar (98%) to Chaco chachalaca 
(Ortalis canicollis). 

The comparative analysis using reference sequences of 
species that were likely to be hunted or consumed at 
the seizure region (Table 1) confirmed these results, as 
revealed by the genetic distance values and the 
neighbor- joining trees. MAM1 presented a genetic dis- 
tance of zero when compared with the only capybara 
cyt-b reference sequence available in GenBank. MAM2 
differed by 0.32% from the five Pampas deer specimens. 
The genetic divergence between the BIRD sample and 
the Chaco chachalaca was 0.9%. These genetic distance 
values were as low as those for the intraspecific genetic 
distances obtained for the species represented for more 
than one specimen reference sequence (average of 1.2% 
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Table 1 Cyt-b reference sequences of species that would likely be hunted or consumed at the seizure region 



Popular name 



Species 



GenBank accession numbers (ID reference sequence) 



Mammal species 
Indian cattle 

European cattle 

Domestic pig 

Spotted paca 

Capybara 

Marsh deer 

Red brocket deer 

Gray brocket deer 

Pampas deer 

Collared peccary 

White-lipped peccary 
Lowland tapir 

Bird species 
Chicken 

Red-throated piping-guan 
Blue-throated piping-guan 
Little tinamou 

Undulated tinamou 
Bare-faced curassow 
Chaco Chachalaca 
Speckled Chachalaca 
Razor-billed curassow 
Rusty-margined guan 
Great tinamou 



Bos indicus 

Bos tourus 

Sus scrofa 

Cuniculus paca 

Hydrochoerus hydrochaeris 
Blastocerus dichotomus 
Mazama americana 

Mazama gouazoubira 

Ozotoceros bezoarticus 

Pecan tajacu 

Tayassu pecan 
Tapirus terrestris 

Gallus gallus domesticus 

Aburria cujubi 
Aburria cumanensis 
Crypturellus soui 

Crypturellus undulatus 
Crax fasciolata 

Ortalis canicollis/Ortallis pantanalensis 
Ortalis gutatta 
Pauxi tuberosa 
Penelope superciliaris 
Tinamus major 



AY126697 (Bindl), EF693799.1 (Bind2), EF061238.1 (Bind3), EF061239.1 (Bind3), 
EF061 242.1 (Bind4) 

FJ971088 (Btau1),GU249572.1 (Btau2), DQ186288.1 (Btau3), DQ186284.1 (Btau4), 
AY952963.1 (Btau5) 

GQ351599 (Sscrl), GU1 3581 9.1 (Sscr2), AM492594.1 (Sscr3), AY237529.1 (Sscr4), 
GU 13571 6.1 (Sscr5) 

AY206574.1 (Cpacl), AY206563.1 (Cpac2), AY206572.1 (Cpac3), 
AY206561.1 (Cpac4), AY206560.1 (Cpac5) 

GU136721.1 (Hhydl) 

DQ789176.2 (Bdicl), DQ789173.2 (Bdic2), DQ7891 75.2(Bdic3), DQ789174.2 (Bdic4) 

DQ789230.2 (Mamel), DQ789201.2 (Mame2), DQ789225.2 (Mame3), 
DQ789209.2(Mame4), DQ789224.2 (Mame5) 

DQ789200.2 (Mgoul), DQ789189.2 (Mgou2), DQ789182.2 (Mgou3), 
DQ789203.2 (Mgou4), DQ789184.2 (Mgou5) 

L48404.1 (Obezl), DQ789199.2 (Obez2), DQ789198.2 (Obez3), 
DQ789191.2 (Obez4), DQ789192.2 (Obez5) 

DQ1 79085.1 (Ptaj'1), DQ1 79079.1 (Ptaj2) ; DQ1 79082.1 (Ptaj3), DQ1 79074.1 (Ptaj4), 
DQ1 79065.1 (Ptaj5) 

AY534303.1 (Tped), U66290.1 apec2), AY726775.1 (Tpec3) 

AF056030.1 (Tterl), GQ259949.1 (Tter2), GQ259923.1 (Tter3), GQ259954.1 (Tter4), 
GQ259936.1 (Tter5) 



HQ122606 (Ggall), AF195628.1 (Ggal2), AY029583.1 (Ggal3), AF3541 71 .1 (Ggal4), 
AF028795.1 (Ggal5) 

AY659799.1 (Acuj) 

AY659798.1 (Acum) 



FJ899152 (Csoul), FJ8991 51 .1 (Csou2), FJ899K 
FJ899147.1 (Csou5) 

AY1 39629.1 (Cundl) 

AY659790.1 (Cfasl), AY141 923.1 (Cfas2) 

AF1 65472.1 (Ocan)/AY659783.1 (Opan) 

AY659782.1 (Ogutl) 

AY354484.1 (Ptubl), AF1 65469.1 (Ptub2) 
AY659804 (PsupD 
AF338707.3 (TmajD 



).1 (Csou3), FJ899150.1 (Csou4), 



for mammals and 1.9% for birds), and lower than the cyt 
Z?-genetic distances (<2%) typical of population and 
intraspecific variation observed in mammal and bird 
species [17,18]. The species identifications of the seized 
samples were also supported by clades with 96-100% 
bootstrap values (Figures 2 and 3). Therefore, even with 
the questionable nature of the Genebank sequences, 
these results provide an indication of the species hunted. 
In a most conservative way, considering the family of 
the organisms hunted, the three cases already can be 
considered to be wildlife hunting crimes. 

Contrary to results with cyt-&, species identification 
was not possible using COL The amplification of COI in 



BIRD, which could be solved using primers designed for 
bird species, was not successful. For the mammal spe- 
cies, the BLAST analysis did not match a reference 
sequence based in our threshold (>98%). At the time of 
the manuscript writing, Genbank and BOLD did not 
contain COI reference sequences for capybara or for 
Pampas deer. Although the COI gene is considered the 
DNA barcoding gene [7], reference sequences for several 
species are still unavailable in these genetic databases. 
Researchers have used the cytochrome b gene for species 
identifications [5,19] since it is one of the better repre- 
sented genes in GenBank [19] and has superior ability 
for separating species when compared with COI [20,21]. 
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Figure 2 Neighbor-joining tree of species that could likely have 
been consumed at the seizure region and MAM1 and MAM2 
meat samples, based on the cyt-6 molecular marker. 



The three species killed are widely reported in hunting 
studies [1,22,23]. The capybara is widely distributed in 
South America and is the largest living rodent, weighing 
around 50 kg [23]. Pampas deer, the most endangered 
Neotropical cervid, formerly occupied a range of open 
habitats such as grassland, pampas and savanna 
(cerrado) in Brazil [24]. Its populations are decreasing 
because of habitat conversion for agriculture and cattle 
farming, hunting and attacks by feral dogs [24]. Its 
former range has been reduced to less than 1% [25]. The 
Chaco Chachalaca is a galliform commonly found in the 
Pantanal of Brazil. Cracid species are becoming rare be- 
cause of hunting, and the loss and fragmentation of suit- 
able habitats [22]. 

In this study the cyt-b molecular marker was suitable 
to diagnose the species killed and solve these criminal 
investigations. The suspected poaching in all three cases 
was confirmed with the identification of three wild spe- 
cies, capybara, Pampas deer and Chaco Chachalaca. In 
Brazil, hunting is a criminal offense, and based on the 



BHD 
Ocan 
Ogut 



-Opan 



- PSIf) 

I Acij 



100 L_ Acun 
99 I — Cfasl 



Cfas2 
r-Ptii>1 
100 L ptub2 
72|Ggal1 
lGgaJ2 



100 

67 



Ggal5 
Ggal3 
Ggal4 



-Tmaj 



— CukJ 

-Csou3 
-Csou4 

Csou2 

Csoul 



-Csou5 



Figure 3 Neighbor-joining tree of species that could likely have 
been consumed at the seizure region and BIRD meat sample, 
based on the cyt-6 molecular marker. 
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results of our molecular analysis, the defendants were 
found guilty and punished with fines; they may still be 
sentenced to prison for a period of 6 to 12 months. 

Abbreviations 

COI: cytochrome c oxidase subunit I; cyt-b: cytochrome b; BOLD: Barcode of 
Life Database; mtDNA: mitochondrial DNA; IBAMA: Institute Brasileiro do 
Meio Ambiente e dos Recursos Naturais Renovaveis; MAM1: sample meat of 
the mammal seized in case 1; BIRD: sample meat of the bird species seized 
in case 2; MAM2: sample meat of the mammal seized in case 2; 
PCR: polymerase chain reaction. 
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